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1. Introduction 



Measurements in a gravitational field are most easily interpreted through the use 
of a system of locally inertial coordinates. For an observer following a timclikc 
worldlinc, Fermi- Walker coordinates provide such a system. A Fermi- Walker 
coordinate frame is nonrotating in the sense of Newtonian mechanics and is 
realized physically as a system of gyroscopes [1] [2] . The metric along the path 
is Minkowskian, with first order corrections away from the path that depend 
only on the acceleration of the observer [3J. If the worldline is geodesic, the 
coordinates are commonly referred to as Fermi or Fermi normal coordinates, 
and the metric is Minkowskian to first order near the path with second order 
corrections due to curvature [4]. Applications of these coordinate systems are 
voluminous. They include the study of tidal dynamics, gravitational waves, sta- 
tistical mechanics, and quantum gravity effects [5] 16] \7\ 18] 19] \W\ fTT] . 

Under general conditions, a timclikc path has a neighborhood on which a Fcrmi- 
Walkcr coordinate system can be defined [12] . Some general results are known 
for expansions of the metric in these coordinates. In [13] Li and Ni derived the 
third order expansion of the metric and second order expansion of the equations 
of motion in Fermi- Walker coordinates for general spacetimes, and in [T3] they 
found expansions of the connection coefficients, metric, and geodesic equations 
in Fermi coordinates to third order, fourth order, and third order respectively, 
and gave an iteration scheme for calculation to higher order. Marzlin investi- 
gated weak gravitational fields in [7] and found the expansion of the Minkowski 
metric with small perturbations to infinite order in Fermi- Walker coordinates. 

For particular spacetimes and special timclikc paths, some explicit transforma- 
tion formulas are known. For example, in [5], exact coordinate transformations 
were constructed for specific paths in de Sitter and Godcl spacetimes, but the 
calculations for those examples were possible only because exact solutions for 
certain spacelike geodesies could be obtained in closed form. This is not possi- 
ble in general. In [15] the coordinate transformation mapping the Kerr metric 
written in standard Boyer-Lindquist coordinates to its corresponding form in 
Fermi- Walker coordinates was approximated for a path with fixed space coordi- 
nates and then generalized to circular paths, but the coordinate transformation 
in the reverse direction was not provided. 

To our knowledge, completely general and easily usable transformation formulas 
to and from Fermi- Walker (and therefore Fermi) coordinates do not appear in 
the literature. This paper fills that gap. We calculate an explicit formula for the 
general transformation, and inverse transformation, from a priori coordinates 
to Fermi- Walker coordinates for arbitrary spacetimes, in the form of Taylor ex- 
pansions. The expansions are valid in sufficiently small neighborhoods of any 
timelike path. In one direction, from a priori coordinates to Fermi- Walker co- 
ordinates, the coefficients for the n + 2 order terms involve n-th order partial 
derivatives of connection coefficients along the given timclikc path in Fcrmi- 
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Walker coordinates. Thus, using the results of [13] and [2], the transformation 
law we provide in this direction is immediately available up fourth order for 
Fermi- Walker coordinates and to fifth order for Fermi coordinates (in the case 
that the timelike path is geodesic). The transformation formula for the other 
direction, given by Theorem 3 below, from Fermi- Walker to a given coordinate 
system, is completely self-contained and is exact in the case that the coordinate 
transformation from Fermi- Walker coordinates is real analytic. We elaborate on 
this matter in Section 4. Our methods are more direct than those of [IS], which 
require the solutions of systems of a large number of equations dependent on 
metric coefficients. Our method may in principle be used to calculate coordinate 
expansions and Jacobians to arbitrarily high order and is completely general. 

In Section 2 we introduce notation and define Fermi- Walker coordinates. In Sec- 
tion 3, using the fact that covariant derivatives of coordinate 1-forms are tensors, 
we calculate the Jacobian matrix for the transformation from Fermi- Walker co- 
ordinates to a priori coodinates. Section 4 gives the general transformations 
laws in both directions. In Section 5 we apply the general method to develop 
coordinate transformation formulas for vector fields, and Killing fields, in par- 
ticular. Section 6 illustrates the use of our formulas with examples. Section 7 
gives concluding remarks. 

2. Fermi- Walker coordinates 

Let (M, g) be a four-dimensional Lorentzian C' n manifold. For convenience 
we assume that n > 4, but if necessary the results that follow may be read- 
ily adapted for smaller n. The Levi-Civita connection is denoted by V, and 
throughout we use the sign conventions of Misner, Thorne and Wheeler [3J. A 
timelike path is a smooth map from an open interval on the real line to M, 
whose tangent vector is timclikc. A vector field X is said to be Fermi- Walker 
transported along a timelike path a if X satisfies the Fermi- Walker equations, 
which in coordinate form are given by, 

F s X a = V s X a + ^pX 13 = . (1) 

Here u is the four- velocity along a (i.e., the unit tangent vector), Q a g = a a u$ — 
u a ap, and a a is the four-acceleration. As usual greek indices run over 0, 1, 2, 3 
and lower case latin over 1,2,3. It is well-known and easily verified that 
Fu(u) = 0, and if vector fields X and Y are Fermi- Walker transported along 
a, the scalar product X^Yp is invariant along a. Thus, a tetrad of vectors, 
Fermi- Walker transported along a and orthonormal at one point on er, is nec- 
essarily orthonormal at all points on the path. Moreover, such tetrads may be 
constructed so that one of the orthonormal vectors is the tangent vector u. 
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Let <t(t) denote the parameterization of a by proper time r, and let &o{t), 
e i( r )i e 2(j~), e3(r) be an orthonormal Fermi- Walker transported tetrad along 
(j, with eo = u. We note that a method for constructing Fermi- Walker tetrads 
along a timelike path from their initial values is given in [TO] . The Fermi- Walker 
coordinates 3 relative to the tetrad on a are given by, 



x° ( ex P CT (T)0 J ej(T)) 



x K 



(2) 

(exp a{T) (Xe j (T))=\ K , 



where exponential map, exp p (z/), denotes the evaluation at afhne parameter 1 
of the geodesic starting at the point p in the spacetime, with initial derivative v, 
and it is assumed that the X 3 are sufficiently small so that the exponential maps 
in (JSJ arc defined. From the theory of differential equations, a solution to the 
geodesic equations depends smoothly on its initial data so it follows from Eq. (f5]) 
that Fermi- Walker coordinates are smooth. Moreover, it follows from [TJ] that 
there exists a neighborhood U of a on which the map (a; , x 1 , x 2 , x 3 ) : U — ► M 4 
is well-defined, and it is a diffeomorphism onto the image of U. We refer to such 
a map as a Fermi- Walker coordinate chart (x A , U) for a. By construction, it is 
a nonrotating coordinate system for the observer a [TJ 02 [18] . 

Let {y a } be an arbitrary coordinate system on M defined on an open set con- 
taining a portion (or all) of the timelikc path a. We refer to {y a } as a priori 
coordinates. We assume that the metric tensor is known in the a priori coor- 
dinates, and that connection coefficients may therefore be readily computed in 
these coordinates. Henceforth, we use Greek indices and lower case Latin indices 
exclusively for the a priori coordinates. In addition, we adopt the convention 
that the indices A, B, C, D, E take the values 0, 1, 2, 3, while the indices /, J, K, L 
are restricted to 1,2,3, and we use these upper case Latin indices exclusively 
for Fermi- Walker coordinates. Following this notation, the Fermi- Walker tetrad 
vectors, eo(r), ei(r), e2(r), ^(t), along a are given by, 

where the right side is evaluated at x° = t and x K = 0. The four- by- four 
matrix e^(r) is thus the restriction to a of the Jacobian of the coordinate 
transformation y a = y a (x A ). The inverse of this Jacobian matrix is given by 

FIt A 

«fL- (4) 

Finally, we mention that the non zero connection coefficients in Fermi- Walker 
coordinates, evaluated on a, are given by [3]: 



4 



(5) 



In the case that a is a geodesic so that a = 0, all connection coefficients on a 
vanish. However, partial derivatives of connection coefficients, with respect to 
Fermi- Walker coordinates, on a are in general not zero. In the case of Fermi co- 
ordinates, these derivatives along with expansion of the metric tensor to second 
order in the space variables, were computed in [4j [17] . Higher order derivatives 
and higher order expansions of the metric in Fermi and Fermi- Walker coordi- 
nates are given in [T31 HH [5] . 



3. The Jacobian 



In this section, we begin by calculating Taylor polynomials, centered at a par- 
ticular point on the timelikc path <t(t), for 1-form fields, in the Fermi- Walker 
variables, x°, x 1 , x 2 , x 3 . Without loss of generality, we may expand about the 
point fi(r = 0). The Taylor expansion for a 1-form field V A (x°, x 1 , x 2 , x 3 ) has 
the form, 

V A (x°,x\x 2 ,x*) =V A+ x*^ + r C - B ^s + (6) 

where, on the right side, V A and its derivatives are evaluated at (r, 0, 0, 0), with 
r = 0. Here and in what follows the ellipsis indicates either an infinite sum 
or a finite sum with remainder (in the case that the field V is smooth but not 
analytic) . 

The 0th order terms in Eq. ([6]) may be calculated directly from the Jacobian 
on it, 

V A (r,0,0,0) = e%(T)V a (ir(T)) (7) 

Formulas for the higher order terms in © may be deduced from the fact that 
covariant derivatives and multiple covariant derivatives of tensors are tensors, 



V B V A 



= 4 e a A V p V a 

= e£ e P B e A VjVpVa 



(8) 
(9) 



with analogous third and higher covariant derivative expressions. Thus, from 
Eq. ([8]) it follows immediately that, 

^§(r, 0, 0, 0) = e A V V a + I^e^, (10) 
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Similarly, from Eq. ©, 



-^^b(t, 0, 0, 0) = el 4 e\ V 7 VpV a + T D ABfi V D + T D AB - 



dx c (11) 

AC v D v B + 1 CB v A V D 

Combining Eq. (fTT]) with Eq. (JTUJ) gives, 



r^VnVs + r£ R v. v n 



^^(r, 0, 0, 0) = e? c 4 e\ V 7 V^ Q + r^, c ef,^ 

+ Y%[ele s D V 1 V 5 + T E DC e%V e ] ^ 

+ Fac 4 4VsVp + tS b e a A e 5 D V a V s 

where T BD c = ^g^r^^) and all terms on the right side are evaluated at <j(t). 
Explicit formulas for nth order partial derivatives with respect to Fermi- Walker 
coordinates of the 1-form field {V A } may be similarly obtained in terms of nth 
order covariant derivatives in the a priori coordinate system {y a } and n — 1st 
and lower order derivatives of the connection coefficients with respect to Fcrmi- 
Walkcr coordinates at the point <j(t). Thus, the Taylor coefficients in Eq. ([6]) 
are given by formulas in the a priori coordinates on a. 



Theorem 1. In a neighborhood of a point ct(t) on the timelike path a, the 
Jacobian of the transformation from Fermi- Walker coordinates (x°, x 1 , x 2 , x 3 ) 
to a priori coordinates (y°, y 1 , y 2 , y 3 ) is given by, 



Ja( x ) — q x a ( x ° ' x ,x ' • z;3 ) 

- P a a- t b i v c P a - pP p n r a 1 

— e A + x <. i AB e c e B e A i p v ? 

4- I T C sf 7 P n (-pn r a , ra _p \ (13) 

2 | e C B e A y 7») Pn i ff3 fir) L Pv,l ) 

4- p a r D 4- r D (v E P a — p 1 p v r a \ 

+ e D l ABfi + 1 AB \ l DC e E e C e D l jv) 

_ r D 5 P-pa _ V D C 1 ra I , 

1 AC e D e B l PS 1 CB e A e D l (UJ) f ' 

Proof. Without loss of generality, take r = 0. The Taylor expansion given by 
Eq. © is valid in some neighborhood -B CT (o) of cr(0). Let a point p 6 ^o-(o) 
have a priori coordinates (y° , y 1 , y 2 , y 3 ) and let V be a 1-form field on -Bq-(o) 
whose components relative to {y a } are V a (y° , y 1 , y 2 , y 3 ) . Corresponding to the 
a priori coordinates (y°, y x ,y 2 , y 3 ) there is a unique set of Fermi- Walker coordi- 
nates x , x 2 , x 3 ). Then by virtue of this correspondence, Eq. ([6|) determines 



G 



a map J which transforms the components V a (y° 7 y 1 ,y 2 ,y 3 ) of V at p to the 
Fermi- Walker components Va (x° , x 1 , x 2 , x 3 ) . 

We may in particular apply the map J to each of the following elements of the 
canonical basis of the contangcnt space at p: 



dy° 
dy 1 
dy 2 
dy 3 



(1,0,0,0) 
(0,1,0,0), 
(0,0,1,0), 
(0,0,0,1) 



(14) 



of dy c 



Eq. ([6|) may be used to compute J(dy a ) to find the Fermi- Walker coordinates 

v 

that is, the a-th row, (dy a / dx° , dy a / dx 1 , dy a / dx 2 , dy a / dx 3 ) , of the 

p 

Jacobian matrix of the transformation y a = y a (x°,x 1 ,x 2 ,x 3 ). In matrix form, 

the A-th component, J2(x), at (x°, x 1 , x 2 , x 3 ) of J(dy a ) is given by, 

p 



J%(x) 



(15) 



Now, setting V v = S° (the delta function) in Eqs. ©, ([TU]). ([T2]l. and ©, so 
that \7pV v — — Tp (a(0)), yields Eq. (fT3|) for the Jacobian matrix. ■ 

Using the results of the next section, the Jacobian matrix given (Ti~5|) , and its 
inverse, may be used to transform vector fields Fermi- Walker coordinates 

4. Transformation of coordinates 

In this section, we use the Jacobian (|13p to find coordinate transformations of 
the form y a = y a (x°, x 1 , x 2 , x 3 ) and x A = x A (y°,y 1 ,y 2 ,y 3 ). Let the Taylor 
expansion for y a be given by, 



y {x , x , x , x ) — yQ -t- bj^x -t- Cj^j^x x -t- d asc*^ x x 



„ABC 



(16) 



where yjf = y a (0, 0, 0, 0) = y a (cr(0)). Taking partial derivatives of both sides of 
Eq. P^|) with respect to x A and comparing with Eq. (fT5|) yields the following 
coefficients to third order, 
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h a 



Or a — V c p a — p@ p n V a 
ZC AB — 1 AB e C e B e A 1 0ri 

m a ABC = el e B e\ (W^ + 1*,!^ - T%,„) (17) 



+ e D^AB,C + ^A_B fiDC e °E ~ ^C^D L jri) 

- r D P s p 13 r a - r D p^ <='' r Q 

1 AC e D e B l fid 1 CB e A e D L m»7 

Remark 1. Eq JTBJ) applied to y a (x°, 0,0,0) gives the expansion <x(t) = 
cr(0) + cr'(0)T + icr"(0)r 2 + • • • . Similarly, the expansion for y a (0, x 1 s, 0, 0) 
gives coordinates for points lying on a spacclikc geodesic orthogonal to <t(t) at 
r = 0. 

For the purpose of inverting the series (JT5J) , we employ of the following notation, 

X a = e a A x A = b\x A (18) 

Y a = y a (x°,x 1 ,x 2 ,x 3 )-y% (19) 

2C Pl = 2c AB ep e 7 = V AB e p e y e c - T 0n/ (20) 

3\D% Xn = m a ABC e A e B x eC 

jna i riQ th/j , j^a i i /-t I j^D A B C a 

— _i <5A,jj ' 1 juA 1 <5r; ' 1 i;* 1 Ar; 1 AB.C e <5 e A e r; e D 

T 1 AS L 6 e A e _D i W +1 DC e J e A e jj e Bj 

[r2 c r^ ^ e ^ + rg B r^ e f e £] e £ (21) 

with analogous definitions for higher order coefficients. Note that the coefficients 
C^ 7 , -D^, etc. are symmetric in the subscript indices. Eq. (|16p can now be 
rewritten in the form of a series that is easily invcrtible, 

r=r + cf^x^x 1 + d% Xv x s x x x 11 + ■■■ (22) 

It is readily verified that the Taylor expansion for X a is, 

X a = Y a - C^Y^yy + (2(7^(7^ - Df Xv )Y 5 Y x Y^ + ■■■ (23) 
Thus, from Eqs. (|18|) and (fT9]) we may write, 



+ e A (2C^ - D^)(y s - y s )(y x - y'W y%) 



(24) 



Remark 2. Eq. (|24|) may be easily recast as an expansion about any fixed 
point cr(ro) on <j(t). This is accomplished by redefining e A = e A (ro) and 
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yg = 2/ q O(t )) in Eqs. ^TfJ) , $TE$i - (JUJ) , ([21, and replacing x° by .t° - r 
in Eq. (j24")) . 

For Theorem 2 below, we now make the assumption that the tangent vector 
d/dy° is timclike in U so that y° = t may be selected as a time coordinate in 
the the a priori coordinate system (t, y 1 , y 2 , y 3 ). We assume further that r is an 
increasing function r(t) of t along a so that a may be parameterized by t. This 
excludes causality violations along a. Employing a standard abuse of notation, 
we write a(t) for this parameterization of a (as opposed to a(r(t))). Similarly, 
the Fermi- Walker tetrad e a (r) may be reparameterized by t and we denote that 
parameterization by e a (t) and Eqs. (|20[) and (|21[) are correspondingly modified. 

Theorem 2. With the notation and assumptions of the preceding paragraph, 
the Taylor expansion for the transformation from a priori coordinates to Fermi- 
Walker coordinates in a neighborhood of a is given by, 

x A (t,y\y 2 ,y 3 )=r(t)S A 

+ e£(t)(y k - y k (t)) - e A (t)C%(y= - yi(t))(y k - y k (t)) 

+ eimctfii - n« k w - vmw - vm)(y k - y k (t)) 
+ ••• , 

where as before k = 1,2, 3, and y$ (t) = y a {a{t)). 

Proof. For a given point p £ U with a priori coordinates (t, y , y 2 , y 3 ), Eq. 
(|24p may be revised so as to be an expansion about the point o~(t). This is 
accomplished by redefining = e^(t) and y$ = y a (a(tj) in Eqs. ([T7| - (|2T|) . 
Since yfj = t = y°, Eq. (|2^|) becomes a polynomial in y 1 , y 2 ,y 3 with coefficients 
that depend on t yielding Eq. (|2")li .B 

Remark 3. For the purpose of numerical computations, the expressions for the 
coefficients in Eq. (f25|) . as well as for Eqs. (fT3"|) . (fl6|) . and ([24]l may be simplified 
by using Eqs. (f3"3")) . (|3^|) . and (|55|) which appear below in the proof of Theorem 3. 

Remark 4. Fermi- Walker coordinates along a(t) determine a foliation of a 
neighborhood of a by space slices, each with constant r = x° coordinate. Given 
the coordinates [t.y 1 ,y 2 ,y 3 ) of a point p near a(t), Eq. |25]) may be used to 
locate the space slice containing p by estimating x . 

We next find a series for the inverse transformation to Eq. (|25p . In Theorem 3 
below, (5^, a ' ) represents the coordinate 4-tuple for dy a given by Eqs. ([14]) . The 
parentheses enclosing the index a indicate that 8j should be understood as a 
(0, 1) tensor in the calculations of the coefficients in (|26|) below, rather than as 
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a (1, 1) tensor. 



Theorem 3. For a point p with Fermi- Walker coordinates (r, x 1 , x 2 , x 3 ) suffi- 
ciently close to the timelike path <t(t). the a priori coordinates (y° , y 1 , y 2 , y 3 ) 
of p are given by, 



y a (r, x\x 2 ,x 3 ) = y a (r, 0, 0, 0) + e% (r)x' h 
+ ^(V A12 ^ ) )e^(r)e^(T).T^x^ 

+ ^(V^V^^Je^W^W^W^x*^ + • • • 

+ ^ ( V Mn • • • V M2 *W )e£ (r) • • • e£ (t)z* • • • x J " + ■■■, 

where y a (r, 0, 0, 0) = <r a (r), and where in each term the sum on the //fc's is 
to be carried out before the sum on the Jfc's. If y a is an analytic function of 
(t, x ,x , x 3 ), Eq. (|26|) is an infinite scries. If not, it should be interpreted as a 
Taylor polynomial as indicated above. 

Before proving Theorem 3 we give the following Corollary, which is an imme- 
diate consequence of Theorem 3. However, we also include an independent, 
elementary proof to which wc refer in the proof of Theorem 3. 

Corollary 1. With the same assumptions as in Theorem 3, y a may be expressed 
directly in terms of the a priori connection coefficients as follows, 



y a (r, x\x 2 lX 3 ) = y a (r, 0, 0, 0) + e%(r)x K 

-\T a M {a{r))e^{r)el{r)x J x K 

1 f I ( 2? ) 

+ ^|2r^(a(r))r^(a(r)) - I^Xr)) J 

x e s I (T)e p J (T)e' y K (T)x I x J x K + ■ ■ ■ 

Proof of Corollary 1. Eqs. © may be expressed as the evaluation at s = 1 
of the solution y a (s) of the initial value problem, 



d 2 y a , ^ Q dyP dyi 



ds 2 



F 



= 



01 ds ds 

y a (0) = a a (r) 

£(0)=«*e*(r) 



(28) 



10 



From the initial conditions, the Taylor expansion for y a (s) has the form, 

y a (s) = a a (r) + e a K (r)x K s + a% ^ + ag ^ + • • • (29) 

Similarly, 

r^(y(s)) = r^(a(r)) + T^ s (a(r))(y s (s) - a s (r)) + ■■■ (30) 

In a standard way, substituting Eqs. (|29|) and ([30]) into the geodesic equation 
([28jl. to solve for coefficients yields, 

y Q (s;r,2; 1 ,a; 2 ,2; 3 ) = cr Q (r) + se%(r)x K 

+ ^{2r^(a(r))r^(a(r)) - T^ s (a(r)) 

x e 5 I {T)e l3 J {T)e' y K (T)x I x J x K + ■ ■ ■ 

The result now follows by setting s = 1. ■ 

Remark 5. Using the methods of the preceding proof, it is straightforward to 
compute higher order terms in Eq. (|2"T|) . For a well- written treatment of the 
analogous development of Riemann normal coordinates see [THJ [20] • 

Proof of Theorem 3. The initial value problem (|28p may be reformulated in 
Fermi- Walker coordinates as, 



(r, 0,0,0) (32) 
(O,^ 1 ,^ 2 ,^ 3 ) 

where r is fixed. The solution is a linear function of the affine parameter s given 
by X (s) = (t, sx 1 , sx 2 , sx 3 ), which together with Eq. (f3"2"| yields, 

Tf J x I x J = (33) 

at any point on a and all choices of (x 1 , x 2 , x 3 ). Since Tfj is symmetric in its 
two lower indices, and since Eq. (|3"3"|) holds for all x, it follows that Tfj = 0, a 
fact already noted in the remarks preceding Eq. (J5J). However, differentiating 



d 2 X J 
~lts T 



dX B dX c 



Be- 



ds 



ds 
X(0) 
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the geodesic equation in (|32[) with respect to s and using the linearity of the 
solution yields, 

Tf JiK x I x J x K = (34) 

at any point on a, and differentiating repeatedly yields the analogous higher 
order identities, 

' ./;./../, ./.•'''' (35) 

where the comma indicates that the connection coefficient is differentiated with 
respect to x j3 , ■ ■ ■ ,x Jn , and the result is then evaluated at any point on a. 

For fixed r, the Taylor expansion for y"(r, x , a; 2 , x 3 ) is given by, 

dv a 

y a (T,x\x 2 ,x 3 ) = y Q (r, 0,0,0) + ^^,0,0,0)^ + •■• (36) 

Let represent the coordinate 4-tuple for dy a given by Eqs. (|14p . as described 
above. Then = Jj^Sj = dy a /dx K gives the components of the same 1- 
form in Fermi- Walker coordinates. Thus, Eq. f36|) may be written as, 

y Q (r, x 1 , x 2 , x 3 ) =y a (r, 0, 0, 0) + v£ a) (r, 0, 0, 0)x K 



(n) 



BV 

+ ^r(TA0,0)x«x^ (37 ) 
1 r) 2 V^ a) 

Combining Eqs. (J7J), ([TO]). (jTTJ) with Eqs. and ((MJ) gives the following 

formulas for the first, second, and third order terms in Eq. (|37p . 



dvi a) 



V ( K a) (r, 0, 0, Q)x K = 8$ e£ (r)x K = e a K (r)x K 
'^L{ T , W )x J x K =^^)ef(r)e^(r)x J x K (38) 

|^(r, 0,0,0)3: VV H^^^)ef{r)ef{r)e^{r)x I x J x K 

The analogous formula for the general nth order coefficients may be deduced 
using Eq. (f3"5|) and the higher order analogs to Eq. © as follows, 
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(Vj„ • • • Vj 2 Vj, )x J • • • .v L = (JL-Vj^ ■ ■ ■ j )x J ••• r h 

d d 



v., .-••v./.i;, ).r J 



dx J ndx J «-i Jn ~ 2 '•'"'V' (39) 



9a; J " • • • dx j2 

where all multiple derivatives and covariant derivatives are evaluated at (t, 0, 0, 0). 



Remark 6. A derivation of the transformation formulas in this paper may be 
carried out in the reverse direction. An alternative derivation of Eq. (|16p . and 
thereafter Eq. fig]), is possible by starting with Eq. ([27]) or ([26]). Eq. (fl6]) 
results by computing Taylor expansions in r of the coefficients of Eq. ([2"T]) and 
collecting terms. For that purpose, Eq. ([T]) applied to e^-(r) may be used to 
compute de^(0)/dT. These derivatives may be expressed in terms of connection 
coefficients in Fermi- Walker coordinates. 

5. Vector Fields 

The expansions of Sect. 3 for 1-form fields may be extended more generally 
to tensor fields. If T is a smooth tensor field defined in a neighborhood of the 
timclikc path a in terms of the a priori coordinates {y a }, then we may write, 

dT 

T(t, x\x 2 , x 3 ) = T(r, 0, 0, 0) + x K _(r, 0, 0, 0) 

1 d 2 T (40) 

+ 2^^W^ (T ' '°' 0) + ---' 

and the methods of Sect. 2 may be applied to compute the Taylor coefficients 
for the components of T . For ease of exposition, and to reveal a computational 
short cut, we specialize to the case that T is a vector field with components 
y 3 with respect to the coordinates {y a }- The contravariant versions of Eqs. 
CO),©,©, and the resulting analogs to Eqs. (11)]), HI]), and (fTJ]) yield the fol- 
lowing result. 

Theorem 4. Let a = a(r) be a timelike path parametrized by proper time 
t = x . Let V = V 1 (y° , y 1 , y 2 , y 3 ) be a vector field defined on an open set 
containing a. Then in a sufficiently small neighborhood of ct, the components 
of V in Fermi- Walker coordinates are given by, 
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V A (r,x\x 2 ,x 3 ) = e^ + x K [ey K V /3 V a -ri c e^] 

+ lx J x K [e-}e^e^V V a -ri DtJ efV s 

ri D (e je^V V a -T%e^) (41) 
- rj D efelV 7 V s + r° K e£e s D V S V<*} 
+ 0(3), 

where the coefficients on the right hand side are evaluated at tr(r) and 0(3) 
indicates terms of order 3 and higher in the space variables. 

In the case of Fermi coordinates along a timclikc geodesic, the connection coef- 
ficients vanish, and the result is a substantially simpler formula. 

Corollary 2. With the same assumptions as in Theorem 4 except that a = <j(t) 
is a timelike geodesic, 

V a (t,x\x 2 ,x 3 ) = e A V~< +x K e A e a K V a V~< 

+ \x J x K (e^efV a V^F -V^e°T A K>J ) (42) 
+ 0(3). 

Killing fields play important roles in identifying a concept of energy in the con- 
text of general relativity, studying event horizons, and for other purposes. If V 
is a Killing vector field, then, V Q V pV 1 = R^a av V v ■ Substituting this into Eq. 
(14"2")) yields the following. 

Corollary 3. With the same assumptions as in Theorem 4 except that a = <j(t) 
is a timclike geodesic and V is a Killing vector field, 

V a {t iX \x 2 ,x 3 ) = e A V> +x K e A e a K V a V> 

+ \x J x K {fj^ K e A R\ a X -^^Y A K ^ (43) 

+ 0(3) 

Since V v = e^e^V 7 *, the coefficients of the second order terms, x J x K , may be 
expressed in Fermi- Walker coordinates as \{R A kjd^ D ~ ^kd jV D )- F rom Eq. 
(21) of [13], we may replace D by L in this last expression, thus omitting the 
D = term, and then rewrite these coefficients of second order terms in Eq. 
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l[I511 as, 



+ I (R A klj + R\kj) V L x J x K (44) 

--U A T J T K V L 

— 2 KJL X x v ) 

where in the last line we have used the curvature symmetry, R\^j = R A l[kj]- 
The right side of Eq. ([44]) may be used to replace the coefficient of the second 
order terms in Eq. (j43j) . We note in particular that if the Killing vector V is 
tangent to the geodesic <t(t), then V L = for L = 1, 2, 3 and the second order 
terms in Eq. (|43[) vanish, i.e., 

V A (t, x\x 2 ,x 3 ) = e A V(a(T)) + x K e A e a K V a V^a( T )) + 0(3). (45) 

We include an illustration of this last result in the next section. 

6. Examples 

In this section we illustrate results of previous sections with four examples. In 
Example 1, we apply Eqs. (|25|) and (|45|) to the case of a circular geodesic orbit, 
with a tangent Killing vector, around the central mass M in Schwarzschild space- 
time. We use the results to compare surfaces of simultaneity in Schwarzschild 
and Fermi normal coordinates. In Example 2 we use Eq. (f2"T)) to calculate 
y a (r, x 1 , x 2 , x 3 ) for a timelike trajectory with fixed space coordinates in the 
de Sitter universe, obtaining to 0(4) Chicone and Mashhoon's result in [5]. 
In the third example, we compute coordinates relative to a Zero Angular Mo- 
mentum Observer frame in Kerr spacetime. In Example 4, we compute a co- 
ordinate transformation to Fermi- Walker coordinates for an observer at fixed 
Boyer-Lindquist coordinates in Kerr spacetime. 

Example 1. We take as the a priori coordinates the usual coordinates in 
Schwarzschild spacetime, those of the Schwarzschild observer, y° = t, y 1 — 
r, y 2 = 0, y 3 = </>, in which the metric is given as, 

ds 2 = - (l - —) dt 2 + . + r 2 (d6 2 + sin 2 9dcf> 2 ). (46) 

I 2A1 

r 

In these coordinates, the circular geodesic orbit around the central mass, with 
radius corresponding to radial coordinate ro, is given by, 

o-(t) = (t,r ,ir/2,(3t), (47) 
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where /3 = y/M/r^. We note that r > 3M is necessary for a timelike geodesic 
orbit, and that stable orbits are possible only for ro > 6A/. Thus, in the lan- 
guage of Theorem 2, y^(t) = r , y$(t) = tt/2, and y%(t) = fit. 

It is easily computed that, 



T 



= r{t)= l r A^L t = at. (48) 



For ease of notation, let, 



X = 1™ 



r a - 2M 



v/r (r -3M)' 
I = r \ 



(49) 



M 



r Q - 3M' 



so that e is the energy per unit mass of a test particle on the orbit, and I is 
its angular momentum per unit mass. A version of Eq. ^ for this example is 

1E1, 



e° = (e, 0,0,-0, (50) 
e 1 = 



ly/X sin(a/3t) cos(a[3t) ero sin(a/3i) \ 

>0, t= , (51) 



c 2 



c 



\ r VX 
(0,0,r ,0), (52) 
3 ( — l\/X cos(a(3t) sm(a/3t) er cos(a/3t)\ 



\ r VX VX 

where the ordering of components is given by (t,r,9,<p). Using Eqs. (|20p and 



1G 



I|25p we readily obtain the functions x A (t, r, 8, <ft) to 0(2), 

x° = at- — r (0 - /%) + ••• , (54) 
e 1 

x 1 = e i( r -r o ) + ^(0-/3i) 

+ ^[r; r (r-r o ) 2 +r^(0-7r/2) 2 + r^(0~/3t) 2 ]+..., (55) 

x 2 = r(6 - tt/2) + • • • , (56) 
e 3 

•x 3 = e 3 (r - r ) + ^r(0 - /3t) 

+ y [r; r (r-r o ) 2 +r ee (0- 7 r/2) 2 + r 00 (0-/3t) 2 ]+..., (57) 

where the r^ 7 are the Schwarzschild connection coefficients evaluated at r = ro 
and 9 = tt/2. ft is straightforward to compute higher order terms, but our 
purpose here is merely illustration of method. 

Following Remark 4, we see from Eq. |54|) that 

\x°~at\ w ^-r|^-/%|. (58) 

The left side of Eq. (|58|) vanishes at points on the circular orbit and increases 
with azimuthal deviation from f3t off the orbit and also with increasing radial 
distance from the central mass. To second order, a point with Schwarzschild 
coordinates (t, r, 9, (f>) = (t, r, 9, (3t) lies on the space slice consisting of all points 
(sufficiently close to a) with Fermi time coordinate x° — at. Thus, ignoring 
higher order corrections, we see that the two-dimensional surface consisting of 
fixed time coordinate t and fixed cf> = (3t is simultaneous both for Schwarzshild 
and Fermi observers, i.e., it lies in the intersection of the Fermi and Schwarzshild 
space slices at Schwarzschild time t. As ro — ► oo, x° — > t and the simultaneous 
events in Fermi coordinates become simultaneous for the Schwarzschild observer 
at spacelike infinity, as expected. 

As an application of Eq. |45|) consider the vector field, 

with e and I given by Eqs. (|4"9"| . The right side of Eq. (|B"9"|) gives the components, 
V a , relative to the a priori - in this case, Schwarzschild - coordinates. Since 
V is a linear combination of the Killing vectors d/dt and d/dcf>, with constant 
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coefficients, it is also a Killing vector. Moreover, on the orbit a(t) given by Eq. 
(|47|) , V is given by, 

^(|,0,0,f). (60, 

Thus, V is the four- velocity of the observer following the circular geodesic a(t). 
Using Eq. we easily obtain the components V A in Fermi coordinates as, 

V(t,x\x 2 ,x 3 ) = I 1, -J^x 3 , 0, x 1 I +0(3). (61) 



The observer along ct(t), using Fermi normal coordinates for measurements, 
may define the energy of a nearby test particle following a geodesic path with 
four momentum pa{t, x 1 , x 2 , x 3 ) to be —V a pa- 

Example 2. The de Sitter metric is given by, 

ds 2 = -dt 2 + e 2m [diy 1 ) 2 + d(y 2 ) 2 + d(y 3 ) 2 ] (62) 

where H is a constant. Consider the timelike path a(t) = (t, y 1 , y 2 , y 3 ) with 
fixed a priori space coordinates (y 1 ,^ 2 ,?/ 3 ). In Fermi coordinates, this path 
is parameterized as <j(t) = (r, 0,0,0). Note that from Eq. (f62|) . r = t. An 
orthonormal tetrad along o~(t) is given by, 

e = (1,0,0,0), (63) 

ei = (0,e- fft ,0,0), (64) 

e 2 = (0,0,e- fft ,0), (65) 

e 3 = (0,0,0, e~ m ). (66) 

This tetrad is parallel transported along the observer's geodesic (Chicone and 
Mashhoon [5])- The connection coefficients are 

T\ t =H, T\ i = e %m H, i = 1,2,3, (for y\y 2 , y 3 ). (67) 

The 0(4) contribution in Eqs. and (J27J) is given by 

^ (V^V^M) e v I e^ K elx I x J x K x L = 

(_ p(«) , r" I or^r 1 j_ r' Q ) r A 

- 4r$rv% + 2r ( CT 1r- M rV) e^MV. (68) 
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Using Eq. ([2"T[) and ([55)1 above, we obtain after simple calculation (a large 
number of terms vanish) 



t = y t 



H 



(HR) 2 (HR) 4 
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0(6) 



-(HRf + 0{h) 



(69) 
(70) 



where HR = H\J {x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 . Our results above agree with Eqs. (26) 
and (27) of ref. [5] to fourth order. The method used here produces just as well 
Fermi- Walker coordinates for arbitrary timelike paths. 



Example 3. The Kerr metric in Boyer-Lindquist coordinates is given by 



ds 2 



p 2 A 



dt A + — shr 



o 2 

udt) 2 + ^-dr 2 +p 2 d9 2 , 



where 



p 2 = r 2 + a 2 cos 2 I 



A = 

£ = [r~ + a 
2Mar 



r 2 -2Mr + a 2 , 



- a A sin 



(71) 

(72) 
(73) 
(74) 

(75) 



We consider a Zero Angular Momentum Observer (ZAMO) (or a Locally Non 
Rotating Frame (LNRF) [21]) at some fixed r and 6. The ZAMO's tetrad vectors 
are 



('{) 



f'i 



<'2 



I S IMar \ 

\/4.o,o | , 



0,0,- 



0,0,0, 



S sin^ e 



(76) 
(77) 
(78) 
(79) 



where the ordering of components is given by (t,r,6,<fi). We now go to the 
equatorial plane, 9 = tt/2, and fix the other coordinates t = 0, r = ro, </> = 0. 
Eqs. (|76 p -(|79 [) . evaluated at that spacetime point give the initial condition, or 
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the r = value, of a Fermi- Walker transported tetrad along the circular path 
of a zero angular momentum observer with tangent vector eo- Using Eq. (|2T[) 
with y° =t, y 1 = r, y 2 = 0, y 3 = 4>, and the above tetrad at r = 0, we obtain 

t = ~T t r4> e r 1 etx 1 x 3 + -- - , (80) 

r = , + e ^_i r ; r (eD 2 (x 1 ) 2 -ir^( e «) 2 (x 2 ) 2 



(9 = |+ e e 2 x 2 - T% r e e 2 e r 1 x 1 x 2 + ■■■ , (82) 

(83) 



, 2 



0^3 r<t> „0„r_l„3 



e?x 



3 4 x 0r c 3 c l" 



where at # = 7r/2, r = ro, and A = A(ro), we have 

r r< 4 ZTa ' ^ rr ~~ a ' (84) 



aM (a 2 + 3r 2 ) _ a 2 - r M 

r 2 A ' "■- r A 



'0 'o 



2_ 4 = r 2 (r -2M) -a 2 Af 
ro' 0r r 2 A 



1 0r - — ' 1 <br ~ Ja ■ l 8b J 



Example 4. In this final example we calculate the Fermi- Walker coordinates 
x A (t,r,6,(f)) to 0(2) for a static observer in Kerr spacetime at fixed radial co- 
ordinate ro and azimuthal angle <^ , in the equatorial plane {0 = 7r/2). Wc 
assume that the observer is outside of the ergoregion so that the four ve- 
locity u = (dt/dr, 0,0,0) is timelike. The worldlinc of the observer is then 
cr(i) = (f, ro, 7r/2, (/>o). A Fermi- Walker transported tetrad along this path was 
computed in [16] . Using the notation of Example 3, this tetrad expressed as 
form fields along the observer's worldline may be written, 

e » - (Vx M ™), ,87, 

e 1 = Jo,^cos(7-^o),0,y=^sin(7-^o) J , (88) 



e 2 = (o,-^Lsin( 7 -0 o ),O, J^cos( 7 -0 o ) ) , (89) 



e 3 



/A w T " ' V X 
(0,0,-ro.O), (90) 



where X is given by Eq. (|4"9"| and 



aM , x 

^7v L (91) 
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The corresponding tetrad vectors are given by, 
1 



eo = 



ei 



e2 = 



e 3 



0,0,0 , 



(92) 



2aA/sin(7-^o) VAcos( 7 -0o) n X \ 

7= > , 0, \ — sin 7 - </>o) , 93) 

r y/XA r V A / 



2aMcos 7 -0 o VA sin 7 -^o /X \ 

. , 0, \/ -T- cos 7 - 0o , (94) 

fo V A / 



roVXA 



1 



0,0, ,0 

ro 



(95) 



Using Eqs. (20]) and ([25]) we find x A (t, r, 6, </>) = ^(i, y 1 , y 2 , y 3 ) as, 



x A (t , y 1 , y 2 , y 3 ) = VX^t + e£ (t) (y fc - y k ) 

- \$ic ef(t)e%(t) - e£{t)r*M ~ vtitf ~ Vo) ^ 
+ 0(3), 

where yj = y k (a(t)) is the fcth component of (t, ro, 7r/2, fa) and is independent 
of t for this example (since k = 1,2,3), and the sum over B and C is carried 
out before the sum on j and k in the second term. The connection coefficients, 
r^., in Schwarzschild coordinates appearing in Eq. are given by Eqs. (|84[) 
- (|86[) . Taking into account Eq. ([5]), the nonvanishing Fermi- Walker connection 
coefficients r^ B are given by 



r o _ r o _ r i _ i _ A/x/Acos(7-0o) 
loi-lio-loo-a - r 2 (r -2M) ' ^° 



r o _ r o _ P 2 _„2 _ A/VA sin(7-0 o ) 
02 20 00 " r 2 (r -2M) ' [ ' 

We note that the reverse transformation, from Fermi- Walker to Boyer-Lindquist 
coordinates, may be carried out to arbitrarily high order in a straightforward 
manner using Eq. (|26|) along with Eqs. (l92|) - ([95]) . 



7. Concluding remarks 



Theorem 3 gives the exact transformation formula from Fermi- Walker coordi- 
nates to arbitrary coordinates in general spacetimcs. In particular, Eq. (|26[) 
gives explicit coefficients for expansions in terms of Fermi- Walker coordinates 
to arbitrarily high order. The transformation in the reverse direction is given by 
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Theorem 2, with the Jacobian given by Eq. (fT3"|) . Explicit formulas for vector 
fields were developed in Sect. 5. The examples in Sect. 6 were chosen only to 
illustrate our formulas, which apply quite generally to timclikc paths in general 
spacetimcs. 

Some generalizations of the results of this paper are straightforward. Extensions 
to Riemannian manifolds and to Lorentzian manifolds of arbitrary dimension 
n ^ 2 are easily carried out. Using the methods of this paper, extensions to 
submanifolds beyond timelike paths are also possible. A version of Theorem 3 is 
readily available for the case of Riemann normal coordinates on a Riemannian 
manifold M by applying Theorem 3 to the product manifold of M with the real 
line, with a suitable product metric. 

Rigorous error estimates for the Taylor polynomials of the coordinate trans- 
formations developed in this paper would be useful in some circumstances, as 
would a rigorous lower bound for the radius of a tubular neighborhood of an 
arbitrary timclikc path on which Fermi- Walker coordinates are valid. Different 
techniques are required for the solution of those problems. 
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